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Abstract 

We discuss a general formalism needed for a unified description of the weak and electro- 
magnetic processes in nuclei which is based on the multipole decomposition of the hadronic 
currents. The use of harmonic oscillator single-particle basis, which is commonly employed 
in many-body nuclear calculations, simplifies the relevant expressions. We present analytic 
formulas for the corresponding radial integrals which enter the transition matrix elements of 
one-body operators in various semi-leptonic nuclear processes. As an example, we apply our 
formalism in order to simplify previous formalism giving the nuclear moments for the Ov/3/3 
- decay. Our results for the radial integrals refer to the ls-Od-lp-Of model space. 
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1 Introduction 



In a unified treatment of the semi-leptonic electroweak interactions in nuclei, like lepton-nucleus 
scattering JE 0|, /3-decay modes ||, charged- lepton nuclear capture capture, etc.) M, neutrino- 
nucleus induced reactions 0, ||, ||, exotic semi-leptonic processes in nuclei (double beta-decay 
[0; H El > £* _ — *■ e ~ conversion, etc. |L0|, ITTj), at first, one takes advantage of the well known 



electromagnetic interactions to probe the nucleus and determine accurately the nuclear charge, 
convection current and spin magnetization distributions fl], |2|, 12]. Afterwards, having the 



nuclear structure uncertainties reduced to a minimum, one can predict the cross sections for a 
variety of semi-leptonic weak interaction processes. 

As is well known, the reaction rates between the initial |i) and a final |/) discrete nuclear 
states of any of the aforementioned processes is written in terms of the matrix elements of a 
specific effective Hamiltonians H e ff as 

r^ f ~ \(f\H eff \i)\ 2 

where the H e ff can be assumed to be constructed within the current-current interaction hypothesis 
by first writing down the leptonic and hadronic (J\) currents. 

The effective hadronic current-density operator J\ which enters the H e ff, contains only strong 
isoscalar (T = 0) and isovector (T = 1) pieces and, in general, it has both vector (V) and axial 
vector (A) pieces as 

{ j x] TM T = p(T) { j x] TM T + pip { jlfM T _ (1) 

(T) (T) 

where f3y and f3 A are polar-vector and axial-vector coefficients involving the relevant couplings 
(the purely electromagnetic interactions involve only vector current density, but the weak inter- 
actions involve in addition axial vector current). In neutral-current processes My = for both 
isoscalar and isovector parts, while in charge changing interactions T = 1 with = ±1. One, 
usually, assumes that for a nuclear target (denoted as (A, Z) with Z, A the proton- and mass- 
number respectively) the states \i) and |/) are characterized by definite angular momentum and 
parity (|a) = | J 77 )). In neutral-current reactions (lepton scattering, \i — e conversion, etc.) the final 
state \f) is, in general, an excited state of the nucleus (A, Z) while in charged- current processes 
\f) denotes an excited state of the (A, Z =p 1) nucleus. 

The nuclear calculations of (f\H e ff\i) proceed usually by decomposing the Fourier transform of 
the hadronic current- density operator into irreducible tensors in strong isospin space 0, |3| . Then, 
spherical Bessel functions ji(x) in central or non-central interaction components are obtained. The 
nuclear transition rates are subsequently determined by the matrix elements of these operators. 
The possibility to obtain closed expressions for the matrix elements of the basic one-body operators 
describing the above mentioned reactions has been examined in previous theoretical studies with 
special emphasis on the use of harmonic oscillator wave functions @, §, || |ll|, [l2|, [13 . 



In Refs. |T0[-||T3[ a method was developed which provides explicit compact expressions for 
some basic nuclear form factors (those needed for the /i~ — > e~ conversion in nuclei, scattering of 
cold dark matter candidates off nuclei, etc.). The advantage of this method is that it allows the 
easy calculation of the required reduced matrix elements by separating the geometrical coefficients 
from the kinematical parameters of the reaction. In this way, the transition matrix elements for 
every value of the momentum transfer q can be evaluated immediately. It is worth remarking that 
such a compact formalism provides a very useful insight to those authors who wish to adopt a 
phenomenological approach and fit the nuclear transition strengths. 

In the present work, we extend the method of Ref. [|j|]-|Hl s° as to provide expressions for the 
reduced matrix elements of any of the seven basic single-particle tensor operators describing the 
above-mentioned semi-leptonic processes 0, §, |6|. We especially deal with the operators involving 



1 



the differential operator V (nabla), which have not been studied previously |TT], |12], [Tj|. To 



this end, we first treat the corresponding radial matrix elements which contain derivatives with 
respect to the coordinate r. Finally, we discuss some representative examples as applications of 
our formalism focusing on the improvement of the previous expressions giving the radial part of 
the nuclear moments in the neutrinoless double-beta decay 0-0. 



2 Multipole decomposition of the hadronic currents 

From a nuclear physics point of view, the hadronic current (J M ) which enters the weak and electro- 
magnetic interaction Hamiltonian H e ff is of primary interest. In our convention the four-current 
operator J M is written as 

J,(r) = (p(r),J(r)) (2) 

where p(r) the density and J(r) the three-current operators. The standard multipole expansion 
procedure (see Appendix A) || ||, applied on the matrix elements of J^(r), leads to spherical tensor 
operators which are given in terms of projection functions involving spherical Bessel functions Jl(^) 
and spherical Harmonics Y^(f) or vector spherical Harmonics Yj^'^ J (f ) as: 

M J M {qv) = 5 L jj L (qr)Y^f), (3) 

M^ 1)J (gr) = j L (qr)Y% 1)J (r), (4) 

where 

Y£ 1)J (r) = 1 £(Lmlq\JM)Y£(r)i q , (5) 

m,q 

e q is a unit vector in the direction of the three- momentum transfer q. The magnitude of q (q = |q|) 
is given from the kinematics of the studied process (see e.g. Ref. |10|1 ). 

Using the projections functions Eqs. (f|) and @, in the case of the polar- vector current ( J\) 
we obtain the multipole operators 



iCjcoul 

JM;TMt 



JJM\TM T 



J dvM J M (qv)p(r) TMT i (6) 
J dr (-VM J M (qr)) ■ J(r) TMr , (7) 



1. 



ff M -,TM T = I dr f -V x M#(gr) J " H*)tm t , (8) 



TTm%m t = J drM#(gr) • J(r) TMT . (9) 

which have well-defined parity and are called as: Coulomb, longitudinal, transverse-electric, and 
transverse-magnetic multipoles, respectively The first three operators have parity (— ) J (normal 
parity operators), while the parity of T™ a9 is (— ) J+1 (abnormal parity operator). The analogous 
axial-vector multipoles are 

M 5 JM;TMt = J drM J u {qv)p{xf TMT , (10) 
L 5 JM ;TM T = ijdr \U/M J M {qv)\ ■ 3{t)* TMt , (11) 



Tf M ;TM T = dv[-Vx MjJ(qv) ) ■ 3(t)* TUt , (12) 



,1 

T?m%m t = J dvM J M J (qv) ■ J(vf TMT . (13) 

The first three axial- vector multipoles have parity (— ) J+1 while j 1 J na95 is a normal parity operator. 
For a conserved vector current (CVC) like the electromagnetic, the longitudinal multipoles can be 
written in terms of the Coulomb ones as 

LjM r ,Tm T (q) = ^Mj Mj . TMT (q) (14) 

(go is the time component of the four-momentum transfer, g M = (go, q))- In this case, the number 
of independent operators resulting from the projection procedure is reduced to seven. 

The isospin dependence of the operators (|6D-(|i~3"D includes the operator given by @ 

r 1, T = 0, M T = 

= J TQ = T3> T = 1, M T = (15) 

{ t± = T^(n ± r 2 ), T = 1, M T = ±1 

The exact form of I^ It is determined by the specific semi-leptonic reaction in question (see the 
Introduction). 

The matrix elements of the seven basic operators of Eqs. (|6|)-(|T3D involve isospin dependent 
form factors (see Appendix A). For this reason, we define seven new operators of which the 
spin-space parts are written (in the first quantization) as |2j, || : 

T( M ee M J M (qr) = 5 LJ j L (p)Y M (f), (16) 
T 2 JM = Z J M (qr) = M#-<r, (17) 

Ti M = S'i(gr) = -i{|vx M#(gr)} • a, (18) 



T, 



JM 



S"i(gr) = {-VM^qv)}.^, (19) 



T( M = A J M (qv) = Mif(qr) • -V, (20) 

q 

Ti M = A'i(gr) = -i{±V x M#(gr)} • V, (21) 



T/ M ee Q^(gr) = M^(gr)o- • -V. (22) 

(hereafter we use the shorter notation for the projection functions of Eqs. (f|) and (H) as Mjj/ and 
Y|/, respectively). For the reader's convenience, in Eqs. (fl6l)-(|22|) we adopt the usual notation 
[[J. Using properties of the nabla operator (V), Eqs. (PD, (P|) and ( |2"1"D can be rewritten as: 

Ti M = E'i = [J]" 1 ! - J^M^ 17 + (J + l^M^" 17 } • o-, (23) 



T( M = S"i = [J]' 1 {(J + 1) 1/2 M J + 1J + J^Mfr 1 '} • tr, (24) 
T 6 JM = A'i = [J]- 1 ! - jV 2 Mti u + (J + 1)V J, M£ 1J } ■ ^V, (25) 

(throughout this work we use the common symbol [J] = (2 J + l) 1 ^ 2 ). By glancing at Eqs. (|16|)- 
(p5[), we conclude that, the seven basic single-particle operators of definite parity are built with 
the aid of four operators 

0( M = M M , O j 2 m = Mt/ • a, (26) 
0( M = M L M J ■ iv, 0{ M = M j m <t ■ iv. (27) 

which contain the projection functions M M and Mjj/ or their products with the nucleon-spin er 
and/or the V operator. 



3 The Single-particle matrix elements 

Many quantities of the semi-leptonic electroweak processes are expressed (to a good approxi- 
mation) in terms of single-particle nuclear matrix elements of the one-body operators Tf M ,i = 
1, 2, ...7 ||, H |fl. These matrix elements, by applying the Wigner-Eckart theorem, are written 
as: 

■>JM\a m \ _ (\h-rri! ( 3l J h 



{j imi \Tr\j2m 2 ) = (-r^ 1 ^^ ^ (28) 

In our notation, the single-particle wave functions are labelled (see Appendix B) with the quantum 
numbers (nlsjrrij). The reduced matrix elements of the operators Of, 2 = 1,2, 3, 4, after applying 
the re-coupling relations, can be written in the closed forms shown below. 

1. For the operators Of, with i — 1,2, the reduced matrix elements (ji\\0\ L,S ^ J \\j2) have been 



compactly written as [11, 13 



(ji\\Ol L ' Si)J \\j 2 ) = (h L hWlsinMHp^h), i = 1,2 (29) 
where the symbols (h L Z 2 ) and U[ s are 



(hLl 2 ) = (-)'i * [ Zl ][L][j a ]({j J (30) 



r /1 i 2 ^i 

W/ s = [Ji][J2][J](S + 1) 1 /\S + 2)^ 1/2 1/2 S- . (31) 

{ ji 32 J ) 

In the case of the Fermi-type operator Of, we must put S 1 = in the 9-j symbol of Eq. ( pT|) while 
in the case of the Gamow- Teller operator 0% , S = 1. 



2. The reduced matrix elements of after some manipulation can be cast in the form 



(jt\\M LJ (qv) • -V||j 2 ) = 5>Z(jij 2 ; J^mh^Hp^h), a = ± (32) 
where the coefficients are 

^Oij 2 ; J) = H-) h+L+h+1/2 [hM[J] ( 2k+ 2 Tl ) 1/2 0i LhTi) 

xWe(h,ji, 1/2, 32, h, J)m(L, 1, J, Z 2 , / 1; / 2 T 1), (33) 



1 



with representing the common 6-j symbol 

W 6 (li,ji,l/2,j 2 ,i 2 ,J) = * X / 2 }. (34) 

3. Similarly, for the reduced matrix element of Of we can write 

{j 1 \\M J {qv)a ■ -V||j 2 > = E^(JU 2 ; J)Wi|^(p)KZ 2 ), a = ± (35) 

where 

BtUrf^J) = ±<Wi/2[Mi2](Zi J2j 2 -l 2 )W e (h,ji,l/2,j 2 ,2j 2 -l 2 ,J). (36) 

From Eqs. (p9|) , (|32|) and fl35|) we notice that all basic single-particle reduced matrix elements 
required for our purposes rely on the following three types of radial integrals: 

(n^lTOK^) = J drr 2 R: ih (r)9np)Rn 2 i 2 (r), a = 0, ± (37) 

with 

= ^(p) = Up) (J- p ± 2j2+ 2p 1±X ) • (38) 

The argument p in Eqs. fl38|) is equal to p = qr. We note that the derivatives with respect to p 
appeared in the radial matrix elements come from the application of the gradient formula in the 
matrix elements of Eqs. (|32|) and ([35]) . 



3.1 Expressions for the Radial Integrals 

It is well known that, for single-particle wave functions with arbitrary radial dependence it is not 
easy to perform analytically the integrations over r in Eq. (|37j). These integrals, however, can be 
simplified in the case when harmonic oscillator basis (see Appendix B) is used []2], |13[. Then they 
take the elegant expressions shown below. 

(i) For the operator 9f(p) in Ref. [HJ it was proved that 

Umax 

(nih\3L(p)\n*h) = e~V /2 £ V = (qb/2) 2 (39) 

fi=0 

Umax = (iVl + N 2 - L)/2, 

where N$ = 2rii + l{ represent the harmonic oscillator quanta of the i t h level. The coefficients 
£^(nilin 2 l 2 ) are given by 

efathn-zh) = G— £ £ n!A mi (n 1 / 1 )A m2 (n 2 / 2 )A M (nL), (40) 

m\=tf>ni2=cr 

with 

n = mi + m 2 + {h + h - L)/2. 

Also, G = b 3 N ni i 1 N n2 i 2 /2, where A^ niil is defined in Appendix B. The other symbols of Eq. (f40"D 
are explained in Ref. 



(ii) The formulation of the radial matrix elements which include differential operators Of, 
which constitute the main task of our present paper, proceeds in a similar manner to that of Eq. 
This leads to the expressions 

(A 97 _l_ 1 -I- 1 \ Tlmax 

t p ± 2 P ) ln2h) = e ' y y {L ~ m £ $w v*, ( 41 ) 

where the geometrical coefficients Cfi n ihn 2 l 2 ', L) are given in terms of those of Eq. (fTO|), i.e. 

(L) = _1 f (n 2 + l 2 + 3/2) 1 /2 £ i( ril i in2 Z 2 + i) + n2 /2 £j(niZm 2 - 1Z 2 + 1), < /i < rw* 
^ ^ 2 \ (n 2 + / 2 + 3/2) 1 / 2 ^_(n 1 W 2 + l), M = *W 

(42) 

+ If (n 2 + Z 2 + l/2)V2 e J( ni / l7l2 / 2 _ i) + ( n2 + l)V2 e J( ni / in2 + l/ 2 - 1), < p < n max 

^ { ' 2\(n 2 + ly^Jmhnz + ll 2 - 1), fi = n max 

(43) 

The value of the index n max in the latter two cases of radial integrals is determined by 

n max = (N 1 + N 2 - L + I)/ 2. 

At this point, it should be noted that, the explicit and general formulas of Eqs. (f40D, ( [42]) 
and fl43|) hold for every combination of the levels (nxl\)j\, {n 2 l 2 )j 2 , and give the geometrical 
(momentum independent) coefficients of the polynomials defined in Eqs. ( p9|) and (|4T|). These 
closed expressions have been constructed by inverting properly the multiple summations involved 
in the corresponding matrix elements (see Ref. [14, 15|), so as the final summation is performed 
over the harmonic oscillator quantum number N = 2n + I. 

In order to show the advantages of the above formalism we discuss below some applications. 



4 Applications 

As mentioned before, our formalism is applicable for an infinite single-particle basis and can be 



applied to any semi-leptonic electroweak process which take place in the field of nuclei [0, [15]. 
In a realistic case one is forced to truncate the basis set. Thus, e.g. the description of the nuclear 
moments for the neutrinoless double beta decay 48 Ca — > 48 Ti can be done using a model space 
including the seven orbitals of the major shells 2hu and 3hu || ||. 

In Tables 1-3 we list the coefficients ej^ and (f needed to evaluate the radial matrix elements 
(nili\01(p)\n 2 l 2 ) in the above (ls-Od-lp-Of) model space. For some additional geometrical coeffi- 
cients e^, which determine the transition matrix elements for scattering of dark matter candidates 
off nuclei in the model space including the orbitals 2p-0h and Oi-lg, the reader is referred to Ref. 
0- 

By exploiting the analytic expressions of the Sees. 2 and 3, in the next subsection we simplify 
the formalism of Ref. || 0] constructed for the calculation of the nuclear moments in neutrinoless 
double-beta decay. 

4.1 Nuclear moments for the neutrinoless double-beta decay 

In the description of the nuclear moments for the neutrinoless double-beta decay the fol- 

lowing radial integrals are needed: 

niijHhn^-q) = (mhljLiqry^k) (44) 



f: 



ni^imhn^n'^n'^uj) = J q 2+K v(q; ^TZl^mhn^; qjK^n'^n'^ q)dq (45) 

where v(q;u) are functions describing the energy of the intermediate neutrinos in various gauge 
models (light-neutrino mixing, etc.) 0, 0, @, 0] • The parameter u is related to the excitation energy 
of the intermediate nucleus. Using the compact formalism presented in the previous section the 
above integrals are simplified and written as follows: 

a. For the radial moments 7Z1{nilin 2 l 2 ; q) working as in Eq. (|39"D we find 



nUmkn^q) = &V /2 e" y £ ejy" (46) 



A 4 

At=0 



which implies that the nuclear moments TZ1(nilin 2 l 2 ; q) are simply obtained with the aid of the 
coefficients s^{nilin 2 l 2 ) of Eq. (|40|). Note that for k = Eq. fl46|) reduces to Eq. (j39l). 

b. By replacing the radial integrals in the right-hand side of Eq. (ffH) with Eq. ( |5D| ) and 
manipulating their product we take 

K^nxhn^q^lin^n'^q) = e"V Ll+L2)/2 E ^ ( 4 ?) 

fi=0 

where 
with 

nimax = n 1 +n 2 + (h +h~ L x )/2 , n 2max = n[ + n 2 + (l[ +l 2 - L 2 )/2 
The coefficients c a in Eq. ([470 are obtained from those of 



c l > = 'E&ef* v (48) 

A=0 

In the last summation e^ 1 = 0, for n lmax < v < n max , and e L *_ v = 0, for n 2max < [M - v < n max . 
By inserting Eq. fl4"T|) into Eq. ( ^3|) and putting q = 2y 1 / 2 /b [see Eq. (Bj^)], we have 

ni^nxhn^An'A^) = jt^ E c m J V) (49) 

where 

A = / u+(Li + L 2 + k+1)/2 
The quantities X A (tu) represent the integrals 

T» = / v{y-u)e- 2y y x dy. (50) 
Jo 

which can be computed after providing the exact form of v(y;u) ||. Equation (f49|) is much 
simpler than Eq. (D.l) of Ref. |§, since the latter does not include a simple summation as Eq. 
( f49"D , but a double one over similar integrals to those of Eq. (|H3]) . 



5 Summary and Conclusions 

Our main concern in the present paper was the general formalism needed for calculating the weak 
and electromagnetic processes in nuclei. Starting from the decomposition of the hadronic current 
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into tensor multipole operators, we investigated the possibility of constructing analytic formulas 
for the matrix elements of the principal semi-leptonic operators. We payed special attention on 
the multipole matrix elements taken between single-particle states and restricted ourselves on the 
formulation of their radial part. 

By utilizing harmonic oscillator basis, we have achieved all types of radial moments related 
to the seven basic one-body operators (involved in semi-leptonic reactions), to be written in the 
form of closed analytic expressions, i.e. as products of an exponential times a polynomial of even 
powers in the momentum transfer q. The coefficients of the polynomials found are in general 
simple numbers (for diagonal matrix elements they are always rational numbers) and they can be 
calculated with simple codes. We applied the above compact expressions in order to simplify the 
previous formalism |7], |9| giving the nuclear moments in the interesting neutrinoless double-beta 
decay process. This helped us to see the advantages of the formalism described in Sees. 3 and 4. 
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Appendix A 

1. The identities which are needed in the multipole decomposition procedure are 

expHq.r) = 4vr £ (-z) J M J M (r)Y* M (e 

JMj 



'90 1 



e qx exp(—iq.r) 



XM J M J { 



:2vr) 1/2 E(-0 J ( 2J + 1 ) 1/2 
j>i 

where A = ± and the unit vector e qx has the spherical components 

e go = e z = q/q 
A 1 , .. 

e <?± = ~^~~^2 ^ ey ' 

The components of the rank-one spherical tensor V M are 

d_ 

dz 



-V x M 
Q 



J J i 

M \ 



(51) 
(52) 



V 



\dx dy . 



(53) 
(54) 

(55) 
(56) 



2. For a single free nucleon (using Lorentz covariance, conservation of parity, time-reversal 
invariance, and isospin invariance) the matrix elements of the vector and axial- vector currents, in 
spin-isospin space, are written as 

(k'A'; l/2m t ,\ J M (0) T M T |kA; l/2m t ) = m(k'A') [if^ + F^a^q, + iF^gJ w(kA) 



x (l/2rrv|jf T |l/2m t ) 



k'A'; l/2m^(0)TM T |kA; l/2m t ) 



zWk'A' 



F A ] l5l^ ~ iF K p ^75^ - Frjf ; 75<7 AB/ g v w(kA) 



x (l/2m t /|4 /T |l/2m t ) 



(57) 



(5? 



The plane- wave single-nucleon states are labelled with 3- momenta k(k'), helicities A (A') and 
isospins l/2mt(l/2m t i). The single-nucleon form factors F^ — F^ (qV), with T = 0, 1, and 
X = 1, 2, S, A, P, T (vector (Dirac), vector (Pauli), scalar, axial, pseudoscalar,and tensor) are all 
functions of the momentum transfer qf,. 



Appendix B 

The radial part of the wave functions in a (three-dimensional isotropic) harmonic oscillator 
potential is written as 

R nl = N nl x l e-* 2 / 2 £ l n +1 / 2 (x 2 ), (59) 
where x = r/b, with b is the harmonic oscillator size parameter and N n i the normalization factor 

2n' 

N " = W(b + I + 3/2) - (60) 

r(x) denotes the known gamma function and C 1 ^ 1 ^ 2 represent the Laguerre polynomials defined 
by: 

C V2 (-) = t A m (-0 * m =f t ^( n l l ^y 2 ) * m ■ (61) 



., ml \ n — m 

m=0 m=0 

In Sect. 3, by writing \n(ll/2)jmj), we mean R n ij(r)[Yi(Q r ) ® Xi/ilia, with X1/2 the Pauli 
spinor. We note that for the harmonic oscillator it holds R n ij(r) = R n i(r). The adopted, in the 
present work, sequence of single-particle states is: Osl/2, 0p3/2, Opl/2, ... 
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Table 1: Coefficients which determine the radial nuclear moments {n{t\\j L{Q r ) rK \n2h) (see Eqs. 
( p9[) and fl46|) ) in the ls-Od and lp-Of model space. 



i n 



nih 


- n 2 h 


L 


n = o 




fi = 2 


fi = 3 


Od 


- Od 


1 

3 


7 
6 

3 
5 


14 
15 

2 

15 


2 
15 




Of 


-Of 


1 


3 
2 


9 
5 


18 
35 


4 

105 






3 


33 
35 


44 
105 


4 

105 








5 


26 
105 


4 
105 






lp 


— lp 


1 


5 
6 


19 
15 


13 
15 


2 

15 


Od 


-of 


2 
4 


3 ff 
5 V 2 

11 [% 
15 V 7 


6 ff 
5 V 7 

2 /2" 

15 V 7 


2 /2 
15 V 7 




Od 


— lp 


2 





3 
5 


2 

15 




Of 


- Od 


2 
4 


3 ff 
5 V 2 

11 [2 
15 V 7 


6 ff 
5 V 7 

2 /2" 
15 V 7 


2 /J 
15 V 7 




Of 


— lp 


3 


3 ff, 
5 V 7 


13 /2" 
15 V 7 


2 /? 
15 V 7 




lp 


-Od 


2 


14 
15 


13 
15 


2 

15 




lp 


-of 


3 


9 /2 
5 V 7 


17 /J 
15 V 7 


2 /2" 

15 V 7 





Table 2: Geometrical coefficients (nih^h, L) which determine the radial integrals 
(n^^Of (p)\n 2 l2) of Eq. (|37|). For details see the text. 



1 1 



nih 


- n 2 h 


L 


n = o 




/i = 2 


/i = 3 


Od 


- Od 


1 

3 


i 

2 
1 

15 


4 
15 

2 

15 


2 
15 




Of 


-Of 


1 


5 
6 


1 

15 


26 
105 


4 

105 






3 


9 
35 


16 
105 


4 

105 








5 


2 
105 


4 
105 






lp 


— lp 


1 


1 
6 


7 
15 


7 
15 


2 

15 


Od 


-of 


2 
4 


7 /2 
6 V 7 

1 /2~ 

5 V 7 


4 /2 
15 V 7 

2 /J 
15 V 7 


2 /z 
15 V 7 




Od 


— lp 


2 


2 
5 


1 

5 


2 

15 




Of 


- Od 


2 
4 


i /y 

15 V 2 

1 [2 
15 V 7 


8 [2 
15 V 7 

2 /2" 
15 V 7 


2 /2~ 
15 V 7 




Of 


— lp 


3 


1 p2 
5 V 7 


7 [% 
15 V 7 


2 pi 
15 V 7 




lp 


-Od 


2 


4 
15 


1 

5 


2 

15 




lp 


-of 


3 


1 [2 
15 V 7 


i Fl 

5 V 7 


2 pi 
15 V 7 





Table 3: Geometrical coefficients C^( n i^i ?7 -2^2, which determine the radial integral 
(niZi|fy + (p)|n 2 £2) °f Eq. (P7|) . For details see the text. 



1 o 



